Long-range forces : atmospheric neutrino oscillation at a magnetized detector 
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Among the combinations — L^, — and — any one can be gauged in anomaly free 
way with the standard model gauge group. The masses of these gauge bosons can be so light that 
it can induce long-range forces on the Earth due to the electrons in the Sun. This type of forces 
can be constrained significantly from neutrino oscillation. As the sign of the potential is opposite 
for neutrinos and antineutrinos, a magnetized iron calorimeter detector (ICAL) would be able to 
produce strong constraint on it. We have made conservative studies of these long-range forces 
with atmospheric neutrinos at ICAL considering only the muons of charge current interactions. 
We find stringent bounds on the couplings aefi^er ^ 1.65 x 10~^^ at 3cr CL with an exposure of 
1 Mton-yr if there is no such force. For nonzero input values of the couplings we find that the 
potential V^fi, opposes and Ver helps to discriminate the mass hierarchy. However, both potentials 
help significantly to discriminate the octant of 6123 . The explanation of the anomaly in recent MINOS 
data (the difference of Am|2 for neutrinos and antineutrinos), using long-range force originated from 
the mixing of the gauge boson Z' of — Lr with the standard model gauge boson Z, can be tested 
at ICAL at more than 5(7 CL. We have also discussed how to disentangle this from the solution with 
CPT violation using the seasonal change of the distance between the Earth and the Sun. 



PACS numbers: 14.60.Pq 

I. INTRODUCTION 

The large hadron collider at CERN will probe the 
extensions of standard model above the electroweak 
scale. On the other hand, some of the extensions 
below the electroweak scale can be probed at the 
neutrino oscillation experiments. The extensions be- 
low the electroweak scale introduce massless or nearly 
massless gauge [1, 2] or Higgs bosons [3-6] and they 
couple with matter very feebly and remain invisible. 
These lead to the existence of new kind of forces ei- 
ther i) generating deviations of gravitational law at 
short distances, or ii) predicting low mass particles 
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whose exchange will induce forces at long distances, 
generally violating the equivalence principle [7, 8]. A 
number of experiments have been searching for these 
new forces. The null results provide bounds on par- 
ticle physics models, gravitational physics, and even 
on cosmological models [8-10]. 

The bounds on these couplings to baryon and/or 
lepton number [11] can be obtained from the testing 
of equivalence principle [12] (the free fall acceleration 
is same for all bodies independent of their chemical 
content). In [13], the author has used this idea to 
establish a bound on the strength of an hypotheti- 
cal vectorial leptonic force and obtained the bound 
on the "fine structure" constant a ^ 10^"*^. See [14] 
for a review. A comparable limit also comes from lu- 
nar laser ranging [10] , which measures the differential 
acceleration of the Earth and the Moon towards the 
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Sun. 

However, one can extend the SM with an ad- 
ditional C/(l) gauge symmetry without introducing 
any anomaly for one of the lepton flavor combina- 
tions: Lg — L^, — Lr, and Lg — L^. The masses 
of the gauge bosons can be so light that the induced 
forces may have terrestrial range. Then the electrons 
inside the Sun can induce forces on the Earth sur- 
face depending on the lightness of the gauge boson. 
These forces on the Earth may also be from supernova 
neutrinos, galactic electrons depending on its ranges. 
These couple only to electron (and neutrino) density 
inside a massive object. As a result, the accelera- 
tion experienced by an object depends on its leptonic 
content and mass; and thus violates equivalence prin- 
ciple. 

The long-range (LR) forces can play role in neu- 
trino oscillation. The relatively stronger bounds than 
those from testing of the equivalence principle have 
been obtained from solar, atmospheric and supernova 
data [15-19]^. The long-range forces due to galactic 
electrons, which can affect the galactic rotation curves 
if the range Rlr ^ Rga\ (where Rgai is the distance 
from the galactic center ^ 10 kpc), is also now very 
tightly constrained [16]. Finally, all these results are 
consistent with the existing bounds on violation of 
equivalence principle. 

The sign of this potential is opposite for neutri- 
nos and anti-neutrinos and hence can lead to apparent 
differences in neutrino and anti-neutrino oscillation 
probabilities without introducing CP or CPT viola- 
tion. For instance, the recently found discrepancy in 
the survival probabilities of i^^s and P^s in the MI- 
NOS experiment [21] has been explained using the 
mixing of Z' boson of — Lr symmetry with the Z 
boson of the SM model; and the required associated 
parameters to explain this anomaly can be found in 



[22]. 

We have studied the long-range forces with atmo- 
spheric neutrinos at the magnetized iron calorimeter 
detector (ICAL) proposed at the India-based neutrino 
Observatory (INO) [23], which can directly measure 
the potential detecting separately j/^s and D^s. In this 
article we have focussed on: 

- the bounds on the couplings of long-range 
forces, 

- how significantly long-range forces modify the 
sensitivity of the measurement of oscillation param- 
eters; particularly, the octant of 2-3 mixing (sign of 
joct _ _ ^go-j ^YiQ mass hierarchy (sign of 
Ami,), 

- test of the explanation of the MINOS data with 
long-range forces [22]. 

- possibility to disentangle the effect of CPT and 
long-range forces ^. 

We have also studied the changes of sensitivity 
to mass hierarchy and to 023-octant with true (input) 
023 values considering no such forces (This was not 
studied in earlier works.). 

The paper is organized as follows. We discuss the 
full three flavor neutrino oscillation in matter with 
long-range potential in sec. II, the details of the anal- 
ysis method in sec. Ill, and the bounds on the cou- 
plings in Sec. IV. The effects of the long-range po- 
tential on discrimination of mass hierarchy and on 
6'23-octant is discussed in sec. V. Finally, the test- 
ing of the explanation of the anomaly in MINOS data 
using long-range force is described in Sec. VI. The 
discussion and conclusion arc given in VII. 



^ The advantages of atmospheric neutrinos to discriminate 
CPT violation from CP violation and nonstandard inter- 
actions have been discussed in [24], 



There are plenty relic neutrinos and antineutrinos in the uni- 
verse (~ 50 /cm-^), which may screen the leptonic charges of 
the of celestial bodies. However, it has been analyzed in 
detail and shown that the screening is impossible [20] 
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II. OSCILLATION IN PRESENCE OF 
LONG-RANGE POTENTIAL 

The electrons inside the Sun generate a potential 
at the Earth by [15]: 



^ES 



1.3xlO-^^eV {ae^,er/lO 



-50\ 



(1) 

where, ae^,er = 9%,eT/^'^'^ Se^.er are the gauge cou- 
plings of Le — Lfj_ and Le — Lt symmetries. k, 10^^ 
is the number of electrons inside the Sun [25] and 
Res ~ 7.6 x 10^^ GcV~^, the distance between the 
Earth and the Sun. 

In a three neutrino framework the neutrino fla- 
vor states |i^q), a — e^^l^t can be expressed as linear 
superpositions of the neutrino mass eigenstates \vi), 
i = 1, 2, 3 with masses : 



Uai\Vi) 



(2) 



U is the 3x3 unitary matrix. The time evolution of 
the flavor states is 
. d 



(3) 



where, [i^q] is the vector of flavor eigenstates, and 

Wof = [We), Wt,), Wr)]- 

The evolution equation in the presence of matter 
and long-range potential is 



where, 
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The matter term A = 2\/2GFneE = 7.63 x 
10-5 eV^ /o(gm/cc) £;(GeV) eV^ Here, Gp, and p 
are the Fermi constant, the electron number density 
and the matter density of the medium, respectively. 
The evolution equation for antineutrinos has the re- 
versed sign for A, Vlr and the phase 5. We have 



numerically solved full three flavor oscillation in pres- 
ence of matter and long-range forces. 

However, to understand the bounds it is easy to 
consider two flavor /z — r oscillation, (which has been 
studied in [15] to constrain the long-range forces us- 
ing atmospheric neutrino data of Super-Kamiokandc 
experiment [26]). To understand the effect of 1-3 mix- 
ing, one needs to consider the changes of the effective 
oscillation parameters due to Vlr. The fJ. — t oscil- 
lation in presence of Ver is governed by the evolution 
equation 




cos 26*23 

sin 26*23 



sin 26*23 



AE 



p- COS 26*9 



14. 



Then, the survival probability of i^^ 



P^^ = l- sin^ 2023 sin^ ^^l^ ^ 



^6) 



(7) 



where, L is the neutrino flight path length. The ef- 
fective mixing angle 023 and Am§2 are related with 
their vacuum quantities by the relations 

S'm226'23 



sin^ 261. 



and 



Alfl23 ~ A™23 



[(^er- cos 2023)2 +Sin2 2023] 
(Cer- cos 2023)' +Sin2 2023)'/' 



(8) 



(9) 

where, ^er = ^^r-- The potential Ver and the corre- 
sponding ^ change sign for 9. Similarly, for Lg ~ 
gauge symmetry the survival probability can be ob- 
tained and they satisfy 

(10) 



III. THE X ANALYSIS 

To evaluate the potential of ICAL with atmo- 
spheric neutrinos we generate events by NUANCE-v3 
[27] and consider only the muon energy and its direc- 
tion (directly measurable quantities) of the events for 
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a conservative estimation. The energy and angular 
resolutions of the muons at ICAL are very high: 4- 
10% for energy and 4-12% for zenith angle, which are 
obtained from GEANT [28] simulation. The ranges 
are due to different energies and different angles with 
respect to the iron layers. These uncertainties are 
very negligible compared to the uncertainties in re- 
constructed neutrinos due to kinematics of the scat- 
tering processes. 

The major uncertainty arises from the particles 
produced in the event other than muon. One might 
expect that consideration of the hadrons for neu- 
trino energy ^ 2 GeV will improve the sensitivi- 
ties substantially. But, in [29], it has been found 
that there is a very marginal improvement on mea- 
surement of 023 and a small improvement on Am|]^: 
<5(Am|i) = 0.02 X IQ-^eV'^. The fact is that the total 
hadron energy is carried out by multiple low energy 
hadrons. The average energy per hadron is 1 GcV 
and the average number of hadrons per event is ^1 2. 
The energy resolution of the hadrons at this energy 
is ^ 80%, and consequently, the neutrino energy res- 
olutions do not improve significantly after adding the 
hadrons ^. 

The is calculated according to the Poisson 
probability distribution with flat uncertainties of the 
oscillation parameters. The term due to the contribu- 
tion of prior information of the oscillation parameters 
measured by other experiments is not added to to 
examine solely the performance of ICAL. The data 
have been binned in cells of equal size in the logj^Q E 
- L'^ * plane, where L ~ 2Rcos9z- The choice of bin- 
ning is motivated by pattern of the oscillation prob- 



ability P(i^f_, 



in the L - £; plane [30]. The 



distance between two consecutive oscillation peaks 
driven by Am^g increases (decreases) as one goes to 
lower L (E) values for a given E (L). The binning 



^ The detection of the neutral current events may also be pos- 
sible, which have no directional information, but the energy 
dependence of the oscillation averaged over all directions can 
contribute to the total in the sensitivity studies. 



of L has been optimized to get better sensitivity to 
the oscillation parameters. To maintain x^ /d.o.f 1 
for Monte Carlo simulation study, number of events 
should be > 4 per cell [31] (as large number of cells 
at high energies have number of events less than 4 
or even zero and they increase the x^ jd.o.f substan- 
tially beyond 1). If the number is less than 4 (which 
happens in the high energy bins), we combine events 
from the nearest cells. 

The migration of the number of events from true 
neutrino energy and zenith angle cells to muon energy 
and zenith angle cells is made using exact energy- 
angle correlated 2-dimensional resolution functions 
[32]. For each set of oscillation parameters, we inte- 
grate the oscillated atmospheric neutrino flux folding 
with the cross section, the exposure time, the target 
mass, the efficiency and the two dimensional energy- 
angle correlated exact resolution functions to obtain 
the predicted data in each cell va L — E plane for the 
analysis. We use the charge current cross section of 
Nuance-v3 [27] and the Honda flux in 3-dimensional 
scheme [33] . The number of bins and resolution func- 
tions have been optimized in [31]. Both theoretical 
(fit values) and experimental (true values) data for 
X^ analysis have been generated in the same way by 
migrating number of events from neutrino to muon 
energy and zenith angle bins using the resolution func- 
tions [34]. 

The systematic uncertainties of the atmospheric 
neutrino fiux are crucial for determination of the os- 
cillation parameters. We have divided them into two 
categories: (i) the overall flux normalization uncer- 
tainties which are independent of the energy and 
zenith angle, and (ii) the spectral tilt uncertainties 
which depend on E and 0^. 

The fiux with uncertainties included can be writ- 
ten as 



^{E,ez)^^^{E) 
X [1 + 5^(1 cos 



E 

1 + 5e logio 

-0.5)]x [l + ,S/„] 



For E <1 GcV we take the energy dependent uncer- 
tainty 5e = 15% and Eq = 1 GcV and for £■ > 10 
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GcV, Se = 5% and Eq = 10 GcV. The overall flux un- 
certainty as a function of zenith angle is parametrized 
by 5z- According to [33] we use 6z = 4%, which 
leads to 2% vertical/horizontal flux uncertainty. We 
take the overall flux normalization uncertainty <5/„ = 
10% and the overall neutrino cross-section uncertainty 
6, = 10%. 

We evaluate the for both normal hierarchy 
(NH) and inverted hierarchy (IH) with vs and Ds sep- 
arately for a given set of oscillation parameters. Then 
we find the total ("^ xl + xl)- 

We have set the inputs of jAm^jl — 2.5 x 
10~'^eV'^, and 5cp = 0. We marginalize over 
Am|2, 6*23, 013 and a. We have chosen the range 
of Amj^ = 2.0 - 3.0 x IQ-^eV^ 623 = 37° - 54°, 
6I13 = 0° - 12.5°, and ae^,er = - 3 x 10"^^. The 
solar parameters are fixed at their best-fit values: 
Am|i = 7.67 x lO-'^eV^ and sin^e^is = 0.312 [35]. 
The effect of Arn^i comes in sub-leading order in the 
oscillation probability for atmospheric neutrinos when 
E ^ GeV and it is very negligible. 

In this analysis we have considered an exposure 
of 1 Mton.year (which is 10 years run of 100 kTon) of 
ICAL and energy range 0.8 - 15 GeV. 



order of magnitude of solar mass squared difference 
than the atmospheric one play here the main role to 
make it significantly stringent. 

The bounds are relatively stronger at the mag- 
netized detector: ae/i.er ^ 1-65 x 10"^'^ at 3a CL. 
For atmospheric neutrinos it comes mainly due to the 
high precision of A7TI32 and this can be understood 
quite well from simple two flavor survival probability 
of and in vacuum with effective Am^j and 623 
for long-range potentials. This becomes possible due 
to the fact that the 1-3 mixing effect is sub-leading (as 
6*13 < 11.38° [36]) and it can be neglected at this mo- 
ment for simplicity. When the potential Vlr comes 
to the play, it tries to change the effective value of 
Am|j (see Eq. 9) and becomes tightly constrained. 
The sign Vlr is opposite for u and D, and ICAL can 
detect them separately. This makes the bounds more 
tighter at ICAL and lessens the difference in bounds 
between and aer, while the difference is substan- 
tially large at non-magnetized detectors. 

V. EFFECTS ON 2-3 SECTOR 
A. Determination of mass hierarchy 



IV. BOUNDS ON COUPLINGS 

For atmospheric neutrinos at the magnetized 
ICAL the bounds on and cter are shown in the 
left pannel in Fig. 1 and its strength to constrain 
them with upper as well as lower bounds are shown 
in the right pannel in Fig. 1 with a typical nonzero 
input (true) value. We have checked that the bounds 
do not change significantly for different combinations 
of true values of oscillation parameters. 

The bounds ^ 5.5 x 10^^^ and ^ 6.4 x 
10~^^ at 90% CL have been obtained from present 
atmospheric neutrino oscillation data [15]. Consider- 
ing both atmospheric and solar data the bounds are 
Qfe^ <, 3.4 X 10-53 a,nd aer ^ 2.5 x lO'^^ ^t 3cr CL 
[16]. The better precision and smaller value by a few 



In Fig. 2 we show the sensitivity to mass hi- 
erarchy for different true values of 623- We show it 
for three cases: assuming no potential for long-range 
forces {a^ji ~ 0) and with a benchmark input (true) 
value for ag^i and aer, respectively. 

For all three cases the sensitivity increases as one 
goes to the higher 023 values and it is a general feature 
for both cases with and without the potential for long- 
range force. This can be understood from the i'^ flux 
at the detector considering the effect of 1-3 mixing. 
Here we assume no potential for long-range force. The 
ratio of the flux at the detector (oscillated) and at 
the source (original) [29]: 

1^ « if(sin2^?23) - 7(^23) (1 - ^^(^13), (11) 

where, K{sm 2623) is an even function of the deviation 
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Input: a=0, IH, 9,3=7.5°, 623=45° Input: a=1.1538-10"^^ IH, 6,3=7.5°, 623=45° 




FIG. 1: The upper bounds of a^^ and a^r for true (input) a = in the left pannel; and the corresponding upper as well as 
lower bounds for input a — 1.1538 x 10^'''' in the right pannel, respectively. 




FIG. 2: The sensitivities to mass hierarchy for different true (input) 623 values for input IH (left) and input NH (right), 
respectively. For each case we have considered true (input) values of olb. = 0, a^^ = 1.1538 x 10^''^, and Oer = 
1.1538 X 10"^^, respectively. 



(symmetric with respect to change of the octant), and 
Pa (013 ) is a function of 6*13 only, and 

/(M^(4-^) (12) 

which increases quickly with 623, so that for r = 3 — 4, 
/(6»23 = 40°) < f{023 > 50°). Therefore, for 
6*23 < 45° the flux Ffj_ has much weaker dependence 
on 013 than for 6*23 > 45°. This is reflected in the sen- 
sitivity to mass hierarchy as the lower limit of ^13 has 
been taken zero during marginalization and no prior 
contribution from other future experiments has been 
considered. However, in future the lower limit will be 
known from other experiments (as very recently T2K 
puts a lower bound [37]); then the 023-dependence will 
be less and the sensitivity to mass hierarchy will also 



be substantially improved [34]. 

The potential Ve/j, opposes, while Ver helps to de- 
termine the hierarchy (see Fig. 2). As the 1-3 mixing 
is small (^13 < 11.2°), the potential Ver decreases the 
effective value of Am|i. Consequently, the resonant 
energy decreases (see Fig. 2 of [32]), where the atmo- 
spheric neutrino flux is relatively large and it helps 
in determination of hierarchy. On the other hand, as 
2-3 mixing is large, Vep, increases the effective value 
of Amg]^ for 023 > 45° (see Eq. 9). In spite of an en- 
hancement in 1-3 mixing (which happens for both po- 
tentials) the final sensitivity decreases for Ve^ due to 
relatively low statistics at the resonant zones. These 
happen for i> with NH and for D with IH. 
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Input: IH, e, 3=7.5° 
Input 623=45° Input 623=40° Input 623=50° 




Input 623=45° 



Input: NH, 6,3=7.5° 
Input 623=40° Input 623=50° 




FIG. 3: The sensitivities to deviation of 2-3 mixing from maximal mixing as well as its octant for both true (input) IH (left) 
and NH (right). For each case we have considered olr. ~ 0, ~ 1.1538 x 10^'''', Qer ~ 1.1538 x 10^'''', respectively. 



B. Determination of octant of 623 



VI. MINOS ANOMALY 



The sensitivity to 2-3 mixing, mainly the devi- 
ation from its maximal mixing {6 — 623 — 45°) and 
the octant (sign of S) are shown in Fig. 3 for three 
cases: i) ulr = 0, ii) ~ 1.1538 x 10~^^, and 
iii) ckeT ~ 1.1538 x 10~^^, respectively. It is impor- 
tant to note here that octant discrimination is signifi- 
cantly improved in presence of long-range potentials. 
This is due to the fact that effective 1-3 mixing is al- 
ways larger for both potentials. For a given hierarchy 
the change of the sensitivity with the change of true 
(input) octant depends mainly on the magnitude of 
(5"^ = 6*23°*^ — 45°, which again strongly depends on 
the potentials (see Eq. 8). The change of the effec- 
tive value of Arrig]^ works here subdominantly (while 
it was a dominating factor for determination of the 
mass hierarchy). Now, for Ver with NH, sensitivity 
to octant determination is better for 623 > 45° than 
023 < 45° as 6°^ is increased due to Ver for neutrinos 
and the flux is two times higher for neutrinos than 
antineutrinos. This is opposite for 14^: octant de- 
termination is better for 623 < 45° than 623 > 45°. 
Similarly, the results with IH can be understood using 
antineutrinos considering the symmetries in Eq. 10. 



The anomaly in recent MINOS data (difference 
in measured values of 2-3 mass splittings and mixing 
angles for v and P [21]), which implies CPT violation 
or signal for non-standard interactions, has been ex- 
plained by long-range potential due to — Lr gauge 
symmetry [22]. The Sun docs not contain any /i or 
r and hence there is no direct bound from neutrino 
oscillation. However, the potential can be induced 
indirectly on the Earth by the mixing of the gauge 
boson Z' of — Lr with the standard model Z bo- 
son [22]: 



3.6 x lO^i^eV 



10- 



-50 



(13) 



The sign changes for P. From the analysis of at- 
mospheric neutrino data at ICAL with an exposure 
of 1 Mton-yr, we find the bound on the coupling 
a^r > 3.2(4.2) X IQ-^i at 3(5)cr CL and the explana- 
tion of MINOS data [22] can be tested at more than 
5cr CL. 



VII. DISCUSSION AND CONCLUSION 



In this paper, we have estimated the conservative 
bounds on long-range forces with atmospheric neutri- 
nos at ICAL: ae^,eT ^ 1-65 x lO'^^ ^t 3cr CL. This 



8 



bounds are significantly stronger than the present 
bounds. 

The CPT violation and the long-range forces, 
which are the candidates for solution of recent 
anomaly in MINOS data, can be discriminated with 
atmospheric neutrinos at ICAL. The distance between 
the Sun and the Earth varies and the difference be- 
tween aphelion and perihelion is about 3%, which 
causes seasonal variation of the long-range forces. 
But, there should not be any such change for CPT 
violation. 

Assuming one by one nonzero input value of these 



couplings it is found that while the potential Ve^ 
opposes; Ver helps to discriminate the mass hierar- 
chy. However, both potentials help to discriminate 
octant of 623- The anomaly of MINOS data, which 
has been explained using long-range force potential, 
can be tested at ICAL with more that 5a CL. 
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